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is  developed  on  the  basis  of  the  density-matrix  method.  The  effect  appears  when  light 
excites  intersubband  electronic  transitions.  LID  manifests  itself  as  the  photocurrent  of  the 
2-d  electron  gas  in  the  well  plane,  which  depends  on  coherent  electron  txmneling  between 
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I.  INTRODUCTION 


The  effect  of  light-induced  drift  (LID)  in  gas  mixtures  suggested  in  Ref.  1  manifests 
itself  as  a  drift  of  absorbing  particles  in  response  to  optical  excitation  with  frequency 
close,  but  not  exactly  equal,  to  the  resonant  frequency  of  an  internal  transition  in  the 
particles.  Characteristically,  the  spectral  contour  of  LID  is  antisymmetric  with  respect 
to  detuning  from  the  internal  transition  frequency.  LID  has  been  observed^  and  studied 
in  detail  later  (see,  e.g..  Ref.  3  and  the  references  cited  therein). 

The  effect  of  LID  is  based  on  two  features:  (i)  velocity-selective  excitation  of  parti¬ 
cles  and  (ii)  dependence  of  the  translational  relaxation  rate  of  the  particles  on  their  in¬ 
ternal  state.  The  velocity-selective  excitation  of  particles  is  produced  due  to  the  Doppler 
shift.  Consider,  for  instance,  optical  excitation  red-shifted  with  respect  to  the  exact  res¬ 
onance.  Then  the  particles  which  move  toward  the  light  source  sense  a  blue  Doppler 
shift  which  compensates  the  radiation  detuning.  Such  particles  are  more  likely  to  un¬ 
dergo  excitation.  As  a  result,  their  collision  rate  with  a  buffer  gas,  and  consequently,  the 
friction  force  acting  upon  them  are  changed.  For  the  particles  moving  in  the  opposite 
direction,  there  is  no  excitation  and  no  change  of  friction.  Compelled  by  the  disbalance 
in  the  friction  forces  for  these  two  groups  of  particles,  the  absorbing  component  moves 
as  a  whole.  The  effect  of  surface  LID  predicted  in  Ref.  4  and  observed  in  5  is  similar  to 
the  original  LID,  but  is  bzised  on  collisions  of  the  absorbing  particles  with  the  walls  of 
the  confining  cell  rather  than  with  the  buffer  gas. 

LID  of  electrons  in  semiconductors  predicted  in  Refs.  6  and  7  and  later  observed* 
relies  on  parallel  electron-energy  bands  in  a  quantizing  magnetic  field,  where  an  electron 
moves  freely  in  the  field  direction,  and  light  excites  transitions  between  the  Landau 
levels.  The  translational  motion  of  electrons  in  the  field  direction  yields  velocity-selective 

C 

excitation  which  can  be  interpreted  in  terms  of  the  Doppler  shift,  as  has  originally  been 
suggested  in  Ref.  9  (see  also  Ref.  10).  Because  the  electron  translationzd  relaxation 
depends  on  which  of  the  Landau  states  the  electron  occupies,  LID  takes  place. 

Another  case  of  parallel  energy  bands  mentioned  in  Ref.  6  is  that  of  electrons  in 
two-dimensional  quantum  structures,  such  as  inversion  layers,  thin  films,  etc.,  where  the 
quantum  state  in  the  perpendiculair  direction  plays  the  role  of  an  internal  state  of  the 


particle.  Provided  the  translation  relaxation  depends  on  this  state,  LID  occurs  along  the 
projection  of  the  light  wave  vector  k  onto  the  plane  of  the  structure.  However,  no  theory 
of  such  an  effect  is  given  in  Ref.  6. 

The  photon  drag  effect  (PDE)  in  semiconductor  quantum  wells  has  been 
suggested^ emd  observed  experimentally.^®  Like  in  Ref.  6,  velocity-selective  excita¬ 
tion  in  Refs.  11  and  12  is  due  to  the  parallel  subbands,  and  a  difference  of  the  lifetimes 
of  the  subband  states  is  invoked,  which  makes  this  type  of  PDE  to  be  similar  to  LID. 
Characteristically,  this  difference  is  due  to  optical  phonon  emission  which  is  energetically 
possible  in  the  excited  subband.* 

The  effect  of  LID  in  single  quantum  wells*^  is  based  on  a  different  structure  of  the 
electron  wave  function  in  the  normal  direction,  along  the  lines  of  Ref.  6.  Specifically,  the 
state-dependent  relaxation  in  Ref.  14  is  due  to  an  increased  intrusion  of  the  electron  wave 
function  into  the  barrier  regions  in  the  excited  subband  (see  also  Ref.  15).  Provided  that 
the  doping  of  the  barrier  and  well  regions  is  different,  the  scattering  rate  is  different  in 
the  ground  and  excited  subbands. 

The  present  paper  is  devoted  to  the  theory  of  the  LID  effect  in  coupled  quan¬ 
tum  wells.  Such  systems  allow  one  to  fully  capitalize  on  the  different  structures  of  the 
wave  functions  in  the  ground  and  excited  states  to  achieve  state-dependent  relaxation, 
which  originates  from  quantum  delocalization  of  the  electron  driven  by  resonant  coherent 
tunneling  between  the  wells.  Because  of  the  importeince  of  polarization  relaxation  (de¬ 
phasing),  the  density-matrix  approached  is  invoked,  similar  to  Ref.  14.  The  qualitative 
description  of  the  effect  and  equations  of  motion  for  the  density  matrix  are  presented 
in  Sec.  II.  Analytical  solutions  and  numerical  examples  are  considered  in  Sec.  III.  The 
results  are  discussed  in  Sec.  IV. 

II.  QUALITATIVE  DESCRIPTION 
AND  DENSITY-MATRIX  EQUATIONS 

We  consider  the  LID  effect  driven  by  electron  delocalization  over  coupled  quantum 
wells.  This  effect  can  exist  in  different  quantum  structures,  where  coherent  tunneling  of 
electrons  between  individual  quantum  wells  is  possible.  For  the  sake  of  definiteness  and 
simplicity,  we  restrict  ourselves  to  the  simplest  of  such  structures,  a  double  qu2mtum  well 
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% 


(Fig.  1).  The  excited  state  |2)  in  the  right  well  (numbered  I  in  Fig.  1)  and,  say,  the 
ground  state  |3)  in  the  left  {II)  well  are  aligned  by  the  bias  field  as  shown,  or  by  the 
corresponding  design  of  the  double  well  to  allow  the  resonant  tunneling  of  electrons.  The 
barrier  regions  are  supposed  to  be  modulation-doped  to  yield  the  electron  population 
of  the  structure  and  allow  pure  electronic  intersubband  transitions,  which  only  will  be 
considered. 

One  of  the  wells,  say  II,  is  supposed  to  be  doped  to  create  neutral  scattering  centers. 
When  an  electron  is  excited,  by  virtue  of  the  resonant  tunneling  it  is  delocalized  over  both 
the  wells.  The  electron  experiences  more  frequent  collisions  in  the  doped  II  well,  resulting 
in  faster  translational  relaxation  in  its  excited  state  than  in  the  overall  ground  state  |1), 
because  state  |1)  is  mainly  localized  in  the  undoped  I  well. 

For  the  system  under  consideration,  LID  can  be  described  ac  follows.  The  hght 
has  frequency  close  but  not  exactly  resonant  to  the  intersubband  transition  from  the  |1) 
state  to  an  excited  state.  Electrons  move  in  the  well  pleine  due  to  their  thermal  motion 
for  the  nondegenerate  electron  gas  or  inside  their  Fermi  surface  for  the  degenerate  gas. 
Those  electrons  for  which  the  Doppler  shift  compensates  the  light  detuning  are  excited, 
delocalize  over  the  two  wells,  and,  consequently,  experience  faster  translational  relaxation, 
i.e.  greater  friction  force  than  the  electrons  moving  in  the  opposite  direction.  Compelled 
by  this  force,  the  electron  gas  moves  as  a  whole. 

We  emphasize  that  the  effect  is  fundamentally  based  on  quantum-mechanical  delo¬ 
calization,  i.e.  on  the  presence  of  an  electron  simultaneously  in  two  wells.  In  this  case, 
the  collisions  of  electrons  in  the  doped  {II)  well  implies  translational  relaxation  of  the 
electron  gas  as  a  whole,  which  is  required  for  LID.  The  queintum-mechanical  delocaliza¬ 
tion  over  two  wells  is  brought  about  by  coherent  tunneling,  the  condition  for  which  is 
that  the  tunneling  amplitude  r  is  much  greater  them  the  polarization-relaxation  rate  F. 
In  the  opposite  cjise,  irj  •C  F,  the  tunneling  is  incoherent  and  the  electron  gas  separates 
into  two  distinguishable  weakly-coupled  components,  one  in  the  /  well  and  the  other  in 
the  II  well.  For  both  the  components,  the  relaxation  is  not  state- dependent  and,  con¬ 
sequently,  LID  is  absent,  albeit  the  concentration  of  electrons  in  the  II  well  may  be  not 
very  small  as  is  the  csise  in  Ref.  16.  In  this  respect,  the  present  effects  resembles  the 
resistance  resonance*^  and,  similar  to  it,  can  also  be  considered  as  an  independent  test 
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of  quantum  mechanics. 


The  present  effect  differs  from  the  phenomena  described  in  Refs.  11  and  14  because 
it  requires  coherent  tunneling  and  quantum  delocalization  of  electrons  between  the  wells. 
This  brings  about  a  different  excitation  spectrum  and  enhanced  sensitivity  to  dephasing 
(polarization  relaxation). 

To  describe  the  effect,  we  have  to  take  into  account  the  processes  of  intersubband 
excitation,  the  population  and  polarization  relaxation  along  with  tunneling  between  the 
wells  and  the  translational  movement  and  relaxation  in  the  plane  of  the  quantum  struc¬ 
ture.  As  an  adequate  technique,  we  will  use  the  quantum  Liouville  equation  for  the 
one-electron  density  matrix  pab  in  the  Wigner  representation.  This  quantity  is  defined 
as  Pab{p,r)  =  Eqexp(iqr)/9„fc(p,q),  where  /)a5(p,q)  =  («J_p_iq«a,p-i-iq)  is  the  Fourier 
transform^®  of  Pa6(p,  r)  over  r,  a  is  the  electron  annihilation  operator,  a  zmd  b  denote  the 
states  in  the  well  (a,  6  =  1,2,3),  r  is  the  coordinate  vector,  and  p,  q  are  the  momentum 
vectors  which  all  lie  in  the  plane  of  the  wells. 

The  Hamiltonicin  of  the  electron  plus  ein  external  optical  field  has  the  form  H  = 
Ea,p(^a  +  «p)«i,p«a,p  +  Ea,6,p,q  ^a6(«l)«l,p+iq«6,p-iq>  ^here  £„  is  the  energy  of  the 
level  in  the  wells,  Cp  is  the  electron  kinetic  energy,  and  Uab{^)  is  the  Fourier  transform 
of  the  perturbation  matrix  element  between  states  |a)  and  |6).  The  independent  nonzero 
matrix  elements  of  the  perturbation  are  U2i{r)  =  — d2i  [Eexp(ikr  —  iQt)  -|-  c.c.],  where  d 
is  the  dipole  operator,  is  the  frequency,  k  is  the  wave  vector,  and  E  is  the  electric-field 
amplitude  of  the  light  wave,  and  U32  =  t,  where  the  tunneling  amplitude  r  coincides  with 
the  transfer  integral  t  of  Ref.  19.  The  independence  of  Cp  from  the  quantum  state  (a) 
implies  that  the  subbands  in  quantum  wells  are  parallel,  which  is  of  principal  importance 
for  the  theory,  as  discussed  in  Ref.  6. 

The  exact  equation  of  motion  for  p,  i.e.  the  quantum  Liouville  equation,  can  be 
obtained  in  a  conventional  way  by  commuting  the  peiir  operator  p-|q  ^ith  H 

and  adding  the  releixation  matrix  ( “collision  integral” )  S  which  arises  from  the  interaction 
of  an  electron  with  the  thermal  bath  represented  by  other  electrons,  scattering  centers 
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and  phonons.  This  equation  has  the  form 

=  [£,p(P,r)l  +  5Zexp(iqr)  (cp+iq  -  Cp.iq)  p(p,q) 

q 

+  5]]  exp(iqr)  [f;(q)p(p  -  iq,r)-p(p+  |q,r)t;'(q)]  +t'5  . 
q 

(1) 

In  the  first  sum  on  the  right-hand  side  of  Eq.  (1),  the  shift  of  p  by  ±iq  takes  into 
account  the  velocity-selective  excitation  needed  for  LID,  and  such  a  shift  describes  the 
recoil  effect  in  the  second  sum.  As  determined  by  I7(q)  in  Eq.  (1),  the  magnitude  of  q 
is  on  the  order  of  the  photon  momentum  k  and  is  much  less  than  the  typical  electron 
momentum  p.  This  justifies  the  lowest-order  expzmsion  in  q  in  Eq.  (1).  The  first  order  in 
q  suffices  for  the  first  sum,  and  neglect  of  the  recoil  effect  is  possible  in  the  second  sum 
in  Eq.  (1),  giving  the  quasiclassical  Liouville  equation 


*  (1^  +  V- Ij;)  p(p,r)  =  [e-f  C/(r),p] -l-iS  , 
where  v  =  de^jdp  is  the  electron  group  velocity. 


(2) 


We  employ  the  strong  collision  model  for  the  population  and  translation  relaxations 
and  the  relaxation-constant  model  for  the  dephasing,  i.e.,  we  adopt  the  following  form  of 
matrix  5; 


5ii  =7/  -  ni)  ,  522  =  7/  (”2°^  -  , 

•^33  =  HI  -  "3^  ,  Sab  =  -^abPab  for  O  ^  6  . 


(3) 


Here  ria  =  Paa  is  the  population  of  the  |a)  state  (a  =  1,2,3),  n^a^  is  the  equilibrium 
population  of  the  |a)  state,  =  [1  +  exp  (fa  +  ^ p  ~  where  ep  is  the  Fermi 

energy,  T  is  temperature*®,  7/  and  7//  are  the  population  relaxation  constants  (collision 
rates)  in  wells  I  aind  //,  and  Fai  is  the  polairization-relaxation  constant  for  the  transition 
|a)  \b),  which  is  expressed  in  terms  of  the  pure  dephasing  rate  ra6  in  the  form  ra6  = 

2  (7/  +  7//)  +  Note  that  the  condition  7/  ^  7//  is  necessary  for  the  existence  of  the 
effect. 


The  structure  of  the  diagonal  matrix  elements  of  5  implies  that  that  an  electron  after 
a  single  collision  (an  event  of  interaction  with  the  thermal  bath)  acquires  the  equilibrium 
populations  It  is  worth  noting  that  collision  integral  5  (3)  does  not  conserve  the 


6 


total  electron  concentration,  which  is  clear  because  Tr(5)  ^  0  for  7/  ^  7//.  This  is  a 
well-known  general  problem  for  the  strong  collision  model.  It  is  possible  to  overcome  this 
problem  and  conserve  the  total  electron  concentration  by  introducing  effective  equilibrium 
concentrations  as  in  Ref.  4.  However,  such  a  modification  does  not  bring  about  any 
substsmtial  changes  of  the  solutions.  There  is  another  possible  interpretation,  namely, 
that  the  total  electron  concentration  may,  indeed,  be  changed  due  to  optical  excitation. 
This  is  a  likely  situation  in  quantum  wells  because  they  are  open  system  in  which  the 
electrons  are  produced  by  outside  donors  and  their  equilibrium  may  be  shifted  by  optical 
excitation.  We  accept  this  interpretation  and  will  not  modify  5  (3). 

To  solve  Eq.  (2),  we  will  use  the  rotating- wave  approximation  (RWA),  which  is 
equivalent  to  retaining  only  terms  in  p  with  frequencies  close  to  the  exciting  light  fre¬ 
quency  n  cind  neglecting  multiple  harmonics.  This  resonant  approximation  is  well  known 
to  be  good  for  the  exciting  wave  not  too  strong,  so  that  the  field  broadening  of  levels  is 
much  less  than  fi.  In  any  case,  if  the  exciting  field  is  strong  enough  to  violate  the  validity 
of  RWA,  then  the  three- level  scheme  of  Fig.  1  is  not  applicable,  and  excitation  to  higher 
levels  and  the  continuum  should  be  taken  into  account. 


In  RWA,  the  space-time  dependence  of  Pah  can  be  isolated  in  the  form  p\2  = 
p\2  expt(kr  — fit),  pi3  =  pi3  expi(kr  — fit),  P23  =  P23.  where  the  amplitudes  p\2,  P23,  and 
Pi3  Jire  slowly  [with  respect  to  exp  z(kr  —  fit)]  varying  functions  of  r  and  t.  The  equations 
for  these  2unplitudes  follow  from  Eqs.  (2)  and  (3): 


^  +  V  -  —  J  m  =  2Im(Gp,2)  +  7/  -ni)  , 

n2  =  -2Im(Gpi2  +  7-^23)  +  7/  -  n2^  , 

ns  =  2Im(rp23)  +  7//  (^3°^  -  ^13)  , 


d 

d  > 

dr  } 

1  P\2  =  iG*  (n2  -  ni)  -}-  irpxs  -  gi2Pi2 

d 

d  > 

dr  } 

\  Pis  =  iG*P23  +  IT* p\2  —  Sri3Pl3  1 

d 

d  \ 

di^"' 

dr ) 

1  p23  =  —IT*  {ns  —  n2)  +  iGpis  —  323^23 

(4) 


Here  gu  =  i  (Q  -  kv  -  €ai) -h  Fja  and  ffai,  =  -iefta  +ra6,  where  a,  6  =  2,3,  £ab  =  €a  - 


is  the  transition  frequency,  and  G  =  (dz)-^2Ez  with  the  ir-axis  is  perpendicular  to  the 
structure  plane  (note  that  jG]  is  the  Rabi  frequency).  We  point  out  that  the  Doppler 
shift  kv  in  these  formulas  is  not  introduced  phenomenologically,  but  rather  appears  as  a 
result  of  a  consistent  theory  derived  in  RWA  from  the  microscopic  Hamiltonian. 


The  strong  collision  model  for  the  collision  integral  5  has  been  successfully  used 
previously  to  describe  LID  in  gas  mixtures^,  surface  LID'*,  and  LID  in  a  single  quaintum 
well**.  However,  there  is  a  problem  in  its  application  to  double  quantum  wells  and, 
generally,  to  systems  with  coherent  quantum  tunneling,  which  we  would  like  to  point  out 
and  discuss  before  using  this  model.  This  problem  exists  even  without  the  light-wave 
field,  i.e.  for  G  =  0,  in  which  case  Eq.  (4)  csm  easily  be  solved  for  the  two  upper  levels 
(a  =  2,3)  giving 

+  2  |rp  r23  (^23  +  Ha)"’  (4”^/"/  + 

1 +  2|r|^r23(£23 +  ^13)  *(777+77*)  (®) 

IP23p  =  (^2  -  ^3)^  kp  (£23  +  r^)”*  • 


The  problem  mentioned  above  is  evident  in  the  case  of  completely  coherent  tunnel¬ 


ing,  i.e.  for  |r|  »  r23.  In  this  case,  states  |2)  and  [3)  are  mixed  to  form  new  stationary 

1  ■ 


states  |4-)  and  |— )  with  energies  £±  =  - 


these  states  n+  and  n_  can  be  found  as  the  eigenvalues  of  p  (5)  in  the  form 


62  +  £3  ±  ^£23  +  ^  kl^) 


.  The  populations  of 


1 

n±=2 


n2  -I-  na  ±  (n2  -  na)  (  1  + 


4kl 


n 

2 


(6) 


'23  ^  23 

The  problem  is  that  in  the  generic  case  the  populations  (6)  do  not  conform  to 
the  equilibrium  because  they  do  not  have  the  required  form  of  the  Fermi  distribution 
n±  =  (1  -f  exp(£i  -f  fp  —  £f’)/r]”k  provided  that  the  populations  do  [see  after 
Eq.  (3)|. 


However,  as  we  show  below,  this  problems  is  eliminated  in  limiting  cases  which  are 
determined  by  the  validity  of  at  least  one  of  the  following  three  inequalities: 


|£23|  »  kl  . 

£21, £31  ^  T,  kl  t 

T»kl  • 
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(7) 

(8) 
(9) 


Let  us  consider  these  three  cases  in  some  detail,  (i)  The  inequality  (7)  means  the 
large  mismatch  of  the  levels  Jind  insures  that  the  mixing  of  states  |2)  auid  |3)  is  small. 
In  this  ceise,  from  Eq.  (6)  we  obtain  n+  «  and  n_  «  i.e.  the  populations 

n±  are  the  equilibrium  ones,  (ii)  From  Eq.  (8)  it  follows  that  both  the  |±)  levels  are 
almost  unpopulated,  conforming  to  the  equilibrium.  (Hi)  In  the  case  of  Eq.  (9),  two 
subcases  are  possible.  First,  if  |£23|  <C  T,  then  from  Eqs.  (5)  and  (6)  it  follows  that 
w  n_  «  as  expected  for  the  equilibrium.  Second,  if  |£23l  >  T,  then  we 

return  to  the  case  of  Eq.  (7). 

The  inequalities  (7)-(9)  do  not  contradict  each  other.  We  point  out  that  any  or  all 
of  these  conditions  may  be  valid  experimentally.  For  instance,  setting  realistic  values  (see 
Ref.  16  for  the  details  of  computing  these  values)  for  GaAs  quantum  wells  of  ~  100  A 
thickness  as  r  0.5  meV,  £21  ~  £31  ~  30  meV,  and  T  ~  5  meV  (60  K),  we  conclude  that 
both  Eqs.  (8)  and  (9)  are  valid.  This  justifies  the  use  of  the  collision  integral  (3). 


III.  SOLUTIONS  AND  NUMERICAL  ILLUSTRATIONS 
From  symmetry  considerations,  the  LID  current  is  directed  along  the  projection  of 
k  onto  the  quantum  structure  pleine  (say,  edong  the  x-direction).  Its  density  j  can  be 
expressed  in  terms  of  the  solution  of  Eq.  (4)  in  two  equivalent  forms: 

j  =  t  j  (m  +  n2  -f  n3)dp/(27r)^  =  e(l  -  jiifn/i)  J  v^n3dp/{2Kf  ,  (10) 

where  Vi  —  dfp/dpz  is  the  x-component  of  electron  velocity  and  e  is  the  electron  charge. 
The  second  expression  in  Eq.  (10)  corroborates  the  fact  that  the  driving  force  of  LID  in 
the  system  considered  originates  from  the  difference  of  the  relaxation  rates  in  wells  I  and 
II. 


The  part  of  the  homogeneous  stationary  solution  of  Eq.  (4)  needed  to  find  the 
current  (10)  can  be  exactly  found  cis 

-I-  (B  -  A)  +  n\^\iA  +  {BC  -  A^) 

D  +  7/5  +  (2  +  7///7/)  {BC  -  A^)  ’ 

(11) 

where  D  =  7/7//  -f  7//  [5  +  2  (C  —  .4)]  and 

A  =  2|rG|"Re(F->)  ,  5  =  2  |rl"  Re  [(^12^13  +  M")  F"*]  , 

C  =  2  )G|^  Re  1^^^135123  +  F  *  ,  F  =  923  +  ^12  |C|^ 
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Let  us  consider  some  limiting  cases  of  the  solution  (11).  In  the  absence  of  tunneling 
(r  — >  0),  ns  =  and,  consequently,  LID  is  absent  as  expected;  the  effect  is  driven 
by  quantum  delocalization  due  to  tunneling.  In  the  optically  nonsaturated  cjise,  which 
is  determined  by  the  condition  \G^  ^  li  we  see  from  Eqs.  (10)  and  (11)  that 

j  oc  |G|^,  i.e.  the  LID  current  is  proportional  to  the  light  intensity,  as  expected.  The 
analytical  expression  for  this  limiting  case  is  cumbersome.  To  simplify  it,  we  set  n^**^  = 
[see  Eq.  (9)  and  the  discussion  after  it]  and  get  from  Eq.  (11) 


(0) 

nz-n\' 


2cr23  +  07/  (^32  +  I'32) 


2|Gr|^ 

a  =  Re(/~')  ,  c  =  Re(gi3523/  *)  ,  /  =  523  ^fl'i25'i3  +  . 


7/7//(£i2 +  n2)  +  2kl  r23(7/+7//) 


(12) 


For  the  large  level  mismatch,  |£32|  r32,  the  LID  current  [see  Eqs.  (11)  and  (12)]  tends 

to  zero,  as  expected  tor  an  effect  driven  by  resonance  tunneling.  It  can  be  easily  seen  from 
Eq.  (12)  that  when  the  temperature  is  relatively  high,  T  e2i>  the  LID  current  vanishes 
by  virtue  of  — >  0.  Since  for  the  nondegenerate  gas  LID  also  vanishes  for  T  — >  0 

because  of  u  — »  0,  it  is  clear  that  there  exists  an  optimum  temperature  To  ~  £21. 

It  is  possible  to  obtain  an  analytical  expression  for  the  LID  current  only  in  a  special 
case,  which  we  consider  below.  We  adopt  optically  nonsaturated  conditions.  Also,  we 
assume  the  effective  mass  m*  to  determine  the  electron  dispersion  by  setting  v  =  p/m*, 
and  suppose  the  electron  gas  to  be  nongenerate  using  the  Maxwell-Boltzman  distribution 
ni”\v)  oc  N  exp  {~£a/T  —  v'^lvf),  where  vt  =  (2T/m*)^  is  the  thermzd  velocity  and  N 
is  the  2-d  density  of  electrons.  Additionally,  we  aissume  that  the  tunneling  is  coherent, 
i.e.  |r| »  ri2,ri3,r23. 


In  this  case  the  LID  profile  (12)  consists  of  two  nonoverlapping  contours  j± 
arising  from  the  excitation  of  the  corresponding  components  |±)  of  the  split  dou¬ 
blet.  Each  of  these  contours  centers  about  the  corresponding  transition  frequency 


^±1  = 


£21  +  £31  i  {  ^23  4-4(r 


+4(r(^j 


.  From  Eqs.  (10)  and  (12)  we  obtain 


=  2<r?/e«A-±f5-'  \{d,),/  (7,-/  -  7,-')  Z"'  [1  -  exp(-ej,/r)]  C{?±)cos*  9  ,  (13) 


where  Z  =  exp(— £ai/T)  is  the  partition  function,  I  is  light  intensity,  is  the  normzil 
component  of  the  electron  dipole  operator,  9  is  the  angle  between  the  light-polarization 
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vector  and  the  ^-direction  (the  normal  to  the  well  plane),  and  K±  are  dimensionless 
coefficients  and  C  is  a  function^^  of  the  complex  axgiunents  defined  as 

=  ±4|r|^77^  ^£23  +4|T|^j  [2r23(e±i  -£21) +  £237/]  ^^23  +  2(T(^r23  (77^  +77/ )]  > 

C(0  =  r-e{^exp(-^^)  [1  +eif(iO]}  , 

=  {Q  —  e±i)  /kvt  +  iT±/kvt  ,  r±  =  ±  [Fu  {£±1  —  £21)  +  ri2  (£±i  —  £31)]  ^£32  +  4  |t1^^ 

(14) 

Each  of  the  contours  (14),  j+  and  j_,  has  the  same  antisymmetric  form  with  respect 
to  its  central  frequency  e±i  as  for  LID  in  a  single  welF"*.  However,  the  magnitudes  of 
these  contours  axe  different,  determined  by  the  K±  coefficients.  In  the  case  of  exact  level 
matching,  £23  =  0,  Eq.  (14)  is  greatly  simplified,  giving 

-K’i  =  7/// (7/ +  7//)  1  F±  =  i  (Fn  +  ri2)  .  (15) 

In  this  case,  K±  and  F±  do  not  depend  on  r  and  on  the  |±)-state  index.  Consequently, 
the  LID  current  (13)  does  not  depend  on  the  interwell-barrier  thickness,  which  is  the 
principal  property  of  the  coherent  tunneling  regime. 

Let  us  consider  numerical  illustrations  for  the  structure  shown  in  Fig.  1  in  the 
effective-mass  model  for  the  nondegenerate  electron  gas  adopting  the  following  realistic 
parameters:  m*  =  0.067m  (m  is  the  free-electron  meiss),  £21  =  30  meV,  7/  =  0.1  meV, 

7//  =  0.2  meV,  and  |r|  =  1  meV.  These  correspond^®  to  the  GaAs  wells  with  a  width 
150  A  for  the  I  well,  and  the  Alo.1Gao.9As  central  barrier  with  a  thickness  80  A.  The 
temperature  T  is  set  in  the  calculations  zis  10  meV  (116  K). 

Figure  2  displays  the  spectral  profiles  of  the  LID  current  for  different  polarization- 
relaxation  rates  in  the  optically  nonsaturated  regime  with  exact  level  alignment  (£23  = 

0).  These  profiles  are  calculated  with  Eqs.  (10)  emd  (12)  in  terms  of  the  drift  velocity 
vj  =  j /eN  normalized  to  vt  and  to  \G\^  (meV^)  to  exclude  the  light  intensity. 

For  the  case  Faj  \t\  (weak  polarization  relaxation,  or  nearly  coherent  tunneling), 
the  LID  profile  [see  Fig.  2(a)]  consists  of  two  nonoverlapping  dispersion- type  contours. 

Both  these  contours  are  antisymmetric  with  respect  to  the  corresponding  centers  which 
eire  positioned  at  the  transition  frequencies  of  the  split  doublet  levels,  £±i .  These  contours 
are  quantitatively  described  by  the  terms  sind  j-  of  Eq.  (13).  As  the  calculations  show. 
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in  the  case  under  consideration  the  LID  profile  does  not  significantly  depend  on  r,  i.e. 
on  the  interwell-barrier  thickness,  as  expected  for  coherent  tunneling  [see  the  discussion 
of  Eq.  (15)]. 

For  intermediate  polarization  relaxation,  Tab  ~  |r|.  Fig.  2(b)  demonstrates  that 
the  contours  axe  overlapped,  and  the  antisymmetry  with  respect  to  their  centers  is  lost. 
The  magnitude  of  the  LID  current  is  decreased  20  times,  while  the  increase  of  ra6  is 
only  fivefold.  This  drastic  decrease  of  vj  is  indicative  of  the  separation  of  the  electron 
gas  into  two  distinguishable  gases  localized  in  the  wells  I  and  II,  as  discussed  above. 
A  further  fivefold  increase  of  Faj  [see  Fig.  2(c)]  causes  complete  disappearance  of  the 
two-contour  structure,  which  is  replaced  by  a  single  dispersion- type  antisymmetric  pro¬ 
file  centered  near  the  transition  frequency  £21  in  the  I  well,  as  expected  for  completely 
incoherent  (stepwise)  tunneling.  The  LID  current  is  further  reduced  twentyfold,  man¬ 
ifesting  nearly  complete  loss  of  quantum  coherence  eind  separation  of  the  electron  gas 
into  two  almost  uncoupled  distinguishable  components.  These  findings  directly  confirm 
the  coherent  quantum  nature  of  the  effect.  Note  that  the  intermediate  relaxation  case 
illustrated  in  Fig.  2(b)  is  often  realistic  experimentally. 

The  dependence  of  the  LID  spect^^ll  profiles  on  the  level  mismatch  £23  calculated 
from  Eqs.  (10)  and  (12)  is  presented  in  Fig.  3,  where  the  series  of  curves  is  shown  with 
£21  =  const  and  £31  being  increased.  The  curve  for  the  case  of  perfect  matching  (£23  =  0) 
is  the  same  as  in  Fig.  2(b).  The  change  in  profile  with  an  increase  of  the  mismatc’  is 
rather  nontrivial.  For  moderate  mismatch,  £23  =  2  meV,  the  left  (“red”)  component  of 
the  profile,  which  corresponds  to  the  transition  ]1)  — >  |— ),  is  broadened  and  blue-shifted, 
and  its  amplitude  is  incresised.  The  last  feature  is  counterintuitive,  because  the  effect  is 
driven  by  the  quantum  tunneling  which  is  suppressed  with  the  level-mismatch  increase. 
The  enhancement  of  the  LID  in  this  ceise  is  the  result  of  the  increased  locBdization  of 
the  |  — )-state  wave  function  in  the  I  well.  This  localization  brings  about  better  overlap 
of  the  wave  functions  in  the  ground  state  ]1)  and  the  excited  state  |— ),  which  results  in 
an  increase  in  the  strength  of  the  dipole  transition  jl)  — »  (— ).  Another  reason  for  the 
LID  enhancement  is  the  increase  of  the  electron  lifetime  in  the  |— )  state,  because  of  the 
electron  localization  in  the  undoped  well  I  with  low  collision  rate.  These  two  factors  bring 
about  the  increase  of  the  |— )-state  population  and,  consequently,  LID,  counteracting  the 
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decrease  in  the  tunneling  caused  by  the  intermediate  mismatch. 

The  right  (“blue”)  component  of  the  profile  (Fig.  3),  which  corresponds  to  the 
transition  |1)  —*  |+),  in  all  cases  is  strongly  diminished  with  the  mismatch.  This  is 
because  the  |+)-wave  function  localizes  in  the  II  well,  and  all  the  factors  mentioned 
above  bring  about  a  decrease  of  LID.  When  the  mismatch  becomes  large  (£32  =  8  meV), 
both  the  left  and  right  components  of  the  LID  spectral  profile  decrease  in  magnitude 
due  to  the  diminished  tunneling.  However,  the  right  component  decreases  much  more,  in 
accord  with  the  discussion  above.  In  the  case  of  the  large  mismatch,  the  left  and  right 
components  of  the  spectral  profile  are  centered  at  the  transition  frequencies  in  uncoupled 
wells,  £21  and  £31,  as  expected. 

The  dependence  of  the  drift  velocity  on  (G|  ,  i.e.  on  the  light  intensity,  calculated 
with  Eqs.  (10)  and  (11)  is  shown  in  Fig.  4.  This  dependence  at  small  intensities  is  linear  in 
accord  with  Eqs.  (12)  and  (13),  then  reaches  its  maximum  and  starts  to  decrease  sharply 
due  to  field  broadening  of  the  levels.  With  the  increase  of  dephasing  (cf.  the  dashed  and 
solid  curves),  the  current  reaches  a  lower  maximum  value  at  higher  intensities,  which  is 
due  to  the  separation  of  the  electron  gas  into  two  distinguishable  components,  as  discussed 
above,  and  diminished  velocity  selectivity  of  the  excitation.  For  the  realistic  values  ri2  = 
Fis  =  1  meV,  the  maximum  LID  current  occurs  at  \G^  «  0.1  meV^,  corresponding  to  the 
light  intensity  /  ~  3  kW/cm^,  which  is  easily  reachable  experimentally  [as  an  estimate, 
we  adopted  (dz)i2  ~  maximum  drift  velocity  is  fairly  high,  uj  w  2  •  10~*Ut. 

For  a  typical  2-d  electron  density  N  ~  2  •  10‘®  cm~^  and  well  width  (in  the  y-direction) 
L  =  0.1  cm,  this  yields  the  total  current  per  double  well  of  7  =  eNLvj  «  1  /iA,  which  is 
indeed  easily  detectable. 


IV.  DISCUSSION 

The  theory  of  the  effect  of  light-induced  drift  (LID)  of  electrons  is  presented  in  this 
paper.  Both  qualitative  and  quantitative  properties  of  the  effect  are  discussed  above  in 
Secs.  II  and  III,  in  particular,  in  connection  with  the  numerical  illustrations  (Figs.  2- 
4).  Here  we  do  not  attempt  to  repeat  this  discussion,  but  rather  pursue  two  alternative 
interrelated  goals,  namely,  to  show  the  place  of  this  effect  among  other  electron- transfer 
and  photocurrent  effects  in  quantum  wells,  2Lnd  to  indicate  the  features  of  the  effect  which 
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may  be  useful  for  its  experimental  observation  and  identification. 

LID  in  coupled  wells  requires  the  tunneling  of  electrons  between  wells.  Among  the 
large  number  of  effects  involving  such  tunneling,  most  are  based  on  any  type  of  tunnehng, 
coherent  or  incoherent  (stepwise).  The  transfer  of  an  electron  is  what  is  observed,  usually 
by  recording  photoluminescence  spectra  (see,  e.g..  Refs.  20  and  21).  In  such  experiments, 
the  specific  mechanism  of  tunneling  manifests  itself  only  through  the  timneling  rate  and  its 
dependence  on  the  barrier  thickness.  In  this  connection,  we  point  out  that  even  incoherent 
tunneling  in  some  realistic  cases  can  yield  an  appreciable  electron  transfer.*®  However,  the 
coherence  of  tunneling  is  an  essential  property  for  time-resolved  experiments,  including 
a  recent  one^^  in  which  coherent  electron  oscillations  have  been  observed. 

The  present  effect  manifests  itself  as  the  lateral  stationary  photocurrent  which  de¬ 
pends  on  coherent  tunneling  between  wells.  The  coupling  of  latercil  and  normal  move¬ 
ment  induced  by  the  optical  excitation  makes  LID  unique.  Once  again,  we  emphasize 
that  coherent  tunneling  is  required,  because  collisions  in  one  (doped)  well  should  transfer 
momentum  to  the  electron  gas  as  a  whole.  In  contrast  to  this,  incoherent  tunneling  brings 
about  separate  existence  of  two  distinguishable  electron  gases  in  the  different  wells,  in 
which  case  the  collision  rate  is  not  state-dependent  and  LID  is  absent. 

There  is  another  effect,  namely,  the  resistance  resonance*^,  which  is  also  based  on 
the  coupling  of  lateral  amd  transverse  motion  and  requires  coherent  tunneling.  However, 
the  effect  of  Ref.  17  is  purely  electric,  there  is  no  light,  and  no  intersubbaind  transitions 
and  excited  electronic  states  participate.  Consequently,  different  tunneling  and  relaxation 
processes  contribute.  There  is  another  important  distinction  of  LID:  when  there  is  no 
coherent  tunneling,  LID  is  absent.  By  contrast,  the  tunneling  changes  resistance  typically 
by  20%  only*^.  Therefore,  distinct  from  the  resistance  resonance,  LID  in  the  coupled 
quantum  wells  is  a  background-free  probe  for  coherent  tunneling.  Note  that  different 
rates  of  polarization  relaxation,  i.e.  different  degrees  of  tunneling  coherence,  correspond 
to  essentially  different  forms  of  the  LID  profiles  (see  Fig.  2).  Therefore,  LID  allows  one 
to  obtain  quantitative  information  on  the  tunneling  coherence. 

Now  we  discuss  the  features  of  LID  in  the  coupled  wells  pertinent  to  its  experimental 
observation.  As  estimated  at  the  end  of  Sec.  Ill,  the  maximum  current  is  achieved  under 
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rather  low  light  intensity  (~  1  kW/cm^)  and  is  comparatively  large,  J  w  1  /xA.  We 
emphasize  that  this  is  the  current  induced  in  one  double  well.  Certainly,  such  a  current 
is  easily  detectable.  The  fundamental  property^'*  of  the  LID  current  in  a  single  quantum 
well  of  being  antisymmetric  with  respect  to  the  detiming  from  the  intersubband-transition 
center  is  absent  for  two  coupled  wells  in  the  generic  case.  However,  as  inspection  of 
Eqs.  (10)  and  (11)  shows,  in  the  case  of  exact  level  matching  (£23  =  0),  the  spectral 
profile  of  LID  is  almost  zmtisymmetric  with  respect  to  the  frequency  £21  =  £31,  i.e.  to 
the  center  of  the  split  doublet.  This  property  can  be  easily  traced  in  Fig.  2. 

Another  distinctive  feature  of  the  effect  is  that  it  changes  its  sign  if  well  /,  and  not 
II,  is  doped.  Of  course,  it  is  not  necessary  to  have  a  different  specimen  with  another 
well  doped.  It  is  sufficient  to  change  the  bias  sign  to  have  the  optical  excitation  in  the 
doped  well.  This  is  equivalent  to  switching  wells  /  and  II  (see  Fig.  1),  which  changes  the 
sign  of  the  current.  It  is  advantageous  experimentally  to  have  a  geometrically  symmetric 
double  well  with  one  of  the  individual  wells  doped.  The  double  well  should  be  biased  so 
that  the  the  ground  state  in  one  well  is  aligned  with  the  first  excited  state  in  the  other 
well  achieving  the  scheme  shown  in  Fig.  1.  When  the  bias  reverses  its  sign,  the  current 
also  does. 

We  point  out  the  importance  of  the  alignment  of  the  excited  level  in  one  well  with 
the  ground  state  in  the  other,  as  shown  in  Fig.  1.  Such  a  scheme  excludes  the  counterfield 
electron  transfer*®  which  may  have  interfered  with  LID  if  two  excited  states  were  aligned. 

The  response-relaxation  time  of  LID,  as  follows  from  inspection  of  Eq.  (4),  is  on 
order  of  the  reciprocal  collision  rates  7/*,  ijj  ■  Consequently,  the  realistic  values  of  this 
time  are  on  a  subpicosecond  scale.  Such  a  fast  response  meikes  LID  in  double  quantum 
wells  promising  for  applications. 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  of  the  coupled  potential  wells  I  and  II.  The  overzill  ground  state  |1) 
belongs  in  the  I  well.  The  excited  state  |2)  in  the  I  well  and  the  ground  state  |3)  in 
the  II  well  are  aligned  by  the  bias.  The  optically-excited  intersubband  transition 
is  denoted  by  the  wavy  arrow. 

Fig.  2.  Spectral  profiles  of  the  normalized  LID-drift  velocity  for  the  polarization-relaxation 
constants  ri2  and  Fia  zis  indicated  in  the  figures,  with  r23  =  1.4ri2.  Other  param¬ 
eters  are  given  in  the  text. 

Fig.  3.  LID  spectral  profiles  plotted  for  different  values  of  the  level  mismatch  £32,  as  shown 
in  the  graph,  with  £21  =  30  meV.  The  relaxation  constants  are  ri2  =  ri3  =  1  meV, 
r23  =  1.4  meV. 

Fig.  4.  LID-drift  velocities  plotted  against  |G|^.  The  quantity  |G1^  is  proportional  to  the 
light  intensity,  and  |Gp  =  1  meV^  approximately  corresponds  to  7  =  30  kW/cm^. 
The  constaints  ri2  and  Fis  are  indicated  in  the  figure,  and  r23  =  1.4ri3.  Note  the 
double  logarithmic  scale. 
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Photon  energy  fl  (meV) 
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